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ABSTRACT 

We study magneto spheric structure surrounding rotating magnetized neutron star 
with nonvanishing NUT (Newman-Tamburino-Unti) parameter. For the simplicity of 
calculations Goldreich- Julian charge density is analyzed for the aligned neutron star 
with zero inclination between magnetic field, gravitomagnetic field and rotation axis. 
From the system of Maxwell equations in spacetime of slowly rotating NUT star, 
second-order differential equation for electrostatic potential is derived. Analytical so- 
lution of this equation indicates the general relativistic modification of an accelerating 
electric field and charge density along the open field lines by the gravitomagnetic 
charge. The implication of this effect to the magneto spheric energy loss problem is 
underlined. In the second part of the paper we derive the equations of motion of test 
particles in magnetosphere of slowly rotating NUT star. Then we analyze particle mo- 
tion in the polar cap and show that NUT parameter can significantly change conditions 
for particle acceleration. 
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1. Introduction 



The theoretical study of radio pulsars can be traced back to the work of Goldreich & Ju- 
lian (1969) where the existence of magnetosphere with charge-separated plasma around rotating 
magnetized neutron stars has been suggested. A spinning magnetized neutron star generates huge 
potential differences between different parts of its surface. The cascade generation of electron- 
positron plasma in the polar cap region (Sturrock 1971; Ruderman & Sutherland 1975) means 
that the magnetosphere of a neutron star is filled with plasma, screening the longitudinal electric 
field. This screening results in the corotation of plasma with a neutron star. Such a rotation is not 
possible outside the light cylinder, thus it forms essentially different groups of field lines: closed, 
i.e., those returning the stellar surface, and open, i.e., those crossing the light cylinder and going to 
infinity. As a result, plasma may leave the neutron star along the open field lines and it is generally 
thought that pulsar radio emission is produced in the open field line region well inside the light 
cylinder (the radius at which the corotation speed equals c). 

Study of plasma modes along the field lines is boosted by the pioneering works of Goldreich 
& Julian (1969), Sturrock (1971), Mestel (1971), Ruderman & Sutherland (1975), and Arons & 
Scharlemann (1979). The subsequent achievements and some new ideas are reviewed by numerous 
authors, for example, by Arons (1991), Michel (1991), Mestel (1992), and Muslimov & Harding 
(1997). Although a self-consistent pulsar magnetosphere theory is yet to be developed, the analysis 
of plasma modes in the pulsar magnetosphere based on the above-mentioned papers provides firm 
grounds for the construction of such a model. 

The existence of strong electromagnetic fields is one of the most importa nt features of rot ating 
neutron stars observed as pulsars. It was shown starting pioneering paper of iDeutschl (I1955f) that 
the electric field is induced due to the rotation of highly magnetized star. The ge neral relativisti c 
effect of dragging of inerti al frames is very important in pulsar magnetosphere iBeskinl (|1990|) : 
Muslimov & Tsyganl (| 19901) and a source of additional electric field of general relativistic origin. 



The field-aligned electric field is driven by deviation of space charge from the Goldreich- 
Julian charge density, which is determ i ned by the magnetic field geometry. Therefore as it has been 
noted by several authors (e.g. IBeskinl (119901) . iMuslimov & Tsyganl (119901) . Muslimov & Harding 
(1997); Dyks et al (2001), Mofiz & Ahmedov (2000)) the general relativistic frame dragging effects 
on the field geometry in plasma magnetosphere of rotating neutron stars is a first-order effect, 
which has to be carefully included in a self-consistent model o f pulsar magnetospher i c stru cture 
and associated el ectromagnetic radiation. It was first shown by IMuslimov & Tsyganl (| 19901) and 
independently by IBeskinl dl990|) that general relativistic effects due to frame dragging are crucial 
for formation of the field aligned electric field and particle acceleration in pulsar magnetosphere. 



At present there is no any observational evidence for the existence of gravitomagnetic monopole, 
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that is of exotic space-time, called NUT space (Newman, Unti and Tamburino lNewman. Tamburino & Unti 



(|1963|) ). Therefore it is interesting to study the electromagnetic fields and processes in NUT space 
with the aim to get new tool for studying new important general relativistic effects which are asso- 
ciated with nondiagonal components of the metric tensor and have no Newtonian analogues. 



Muslimov & Tsyganl (119921) initiated the detailed general relativistic derivation of magneto- 



spheric electroma gnetic fields around ro t ating magnetized neutron star. In this paper, we attempt 



to extend work of iMuslimov & Harding! (119971) by including NUT parameter in studying particle 
acceleration along the open field lines of a rotating neutron star. In §[2]general relativistic equations 
describing the electrodynamics of a rotating NUT star are formulated. The equations are rewritten 
in the frame of reference corotating with the rotating NUT star. A detailed analysis of Goldreich- 
Julian charge density is also performed. A general equation governing electrostati c potential in the 



magne tosphere is derived. As it was shown by several authors, in particular by ISakai & Shibata 



(|2003l) . that the general relativistic effect on particle motion in pulsar magnetosphere can not be 
neglected. Due to this reason in § |3] we are solving equations of motion for the charged particles 
for the region near the magnetic pole just above the stellar surface. From study of the particle 
motion along the field lines it has been shown that the effect of NUT parameter, coupling with the 
frame-dragging effect, plays noticeable role in the particle's dynamics. In §@]we conclude our 
findings and discuss them for further investigations. 

Throughout, we use a space-like signature (— , +, +, +) and a system of units in which G = 
1 = c (However, for those expressions with an astrophysical application we have written the speed 
of light explicitly.). Latin indices run 1, 2, 3 and Greek ones from to 3. 



2. Plasma Magnetosphere of Slowly Rotating NUT Magnetized Star 



As it was shown in the paper IMuslimov & Tsyganl (|1992|) . from the system of Maxwell equa 



tions, assuming the magnetic field of a neutron star to be stationary in the corotating frame, the 
following Poisson equation for the scalar potential $ can be derived 

V- (^V$) = -4tt( P - p GJ ) , (1) 

where N = (1 — 2M/V) 1 / 2 is the gravitational lapse function, M is the total mass of the star, J is 
the angular momentum of the star, p — pcj is the effective space charge density being responsible 
for production of unscreened parallel electric field, pcj is Goldreich-Julian charge density to be 
discussed here. 



In a pioneering work, Goldreich & Julian (1969) have shown that a strongly magnetized, 
highly conducting neutron star, rotating about the magnetic axis, would spontaneously build up a 
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charged magneto sphere. The essence of the argument is that it imposes a charge magnetosphere 
that is subject to enormous unbalanced electric forces being parallel to the magnetic field B. Gol- 
dreich & Julian (1969) hypothesized that a far better approximation for the magnetosphere would 
be obtained by shorting out the component of electric field E along B charges originating in the 
star. The magnetospheric charges that maintain E ■ B are themselves subject to the E x B drift that 
sets them into corotation with the star. Here we analyze Goldreich-Julian charge density in general 
relativity in spacetime of slowly rotating NUT star. An expression for the Goldreich-Julian charge 
density p GJ can be calculated in terms of the vector g { = —goi/goo through the formula 

Pgj = -^V(iVg x B) . (2) 



In the metric of the slowly rotating star with the non- vanishing NUT-parameter 



ds 2 = -N 2 dt 2 + N~ 2 dr 2 + r 2 (d9 2 + sin 2 9d(f) 2 ) - 2 (^AN 2 l sin 2 ~ + ur 2 sin 2 9^j d^dt (3) 



vector g looks like following 



g 



N 2 



[CI — ou) x r 



4/iV 2 sin 2 | 
r 2 sin 2 9 



-z x r 



(4) 



This is a stationary axially symmetric solution of vacuum Einstein field equations with three pa- 
ram eters. The metric Q is the lin e ar approximation of the Kerr-Taub-NUT metric (see, for exam- 
ple, |Da^ri£h_&Tjar^uioy| (|2002j), iBini et alj (12003b ) in the specific angular momentum a = J/M 
and the gravitomagnetic monopole I. First term in the right hand side of equation © takes into ac- 
count the well-known effect of dragging of inertial frames of reference (the Lense-Thirring effect) 
with the angular velocity uj = 2aM/r 3 . Here f2 is angular velocity of rotation of star . Obvi ously 
when / = 0, the metric © reduces to the external Hartle-Thorne one JHartle & Thome] Jl968^. On 



other h and when u vanishes the space-time reduces to the NUT one (|Newman. Tamburino & Unti 
(fl963k 



To obtain this equation we remember that in the Cartesian coordinates r = (x, y, z), x 



r sin 9 cos 4>, y = r s i n 9 sin 



r cos 9 and with the z 



place (see lBini et all (12003^ 



fz x r) ■ dr 



2 • 2 

r sm 



(0, 0, 1) the following relation takes 

(5) 



Inserting © into the expression © for the Goldreich-Julian charge density p GJ we obtain 



Pgj 




K 



LI 



e\ 1 4 sin 



2 e 



rj J r] 2 sin 9 



u x B 



(6) 
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where we denoted L = cl/VLR 2 , R is the stellar radius, u = f2xr, rj — r/R is the dimensionless 
radial coordinate, parameter k = ej3, e = 2M/R is the compactness parameter, and (3 = I/I is 
the stellar moment of inertia in the units of 7 = MR 2 . 

In the assumption of the dipole-like configuration for the stellar magnetic field the nonvan- 
ishing components of magnetic field B measured by ZAMO (zero angular momentum observers) 
with four- velocity u a 

u a = {-N 2 , 0,0,0} (7) 
take the form (see Muslimov & Tsygan ( 19921) ) 



B' 



B n ^Mrr]- 3 COS0 , 



7(1) 



where 



-B N 
2 ° 



In 1 



J(rj) 3 

! 7(T) + (i-^)/(i) 



rj 3 sin 9 



s 

2r) 



(8) 



(9) 



and B = 2/i/i? 3 is the Newtonian value of the magnetic field at the pole of star, hat labels the 
ort honormal components, n is the magnetic moment. The solution of this type was first obtained 
by iGinzburg & Ozernoyl (|1964|) and then reproduced by the number of authors. 



The polar angle 9 of the la s t open magnetic line as a function of r\ will look like (see lMuslimov & Tsygan 



(11991b or lMuslimov & Tsvganl (1 19911V ) 



= sin" 1 



1/2 



sin 9r 



e 



sin 



R 



RLcf(l) 



1/2 



(10) 



Parameter 6o has the meaning of the magnetic colatitude of the last open magnetic line at the 
stellar surface, Rlc = 1S me light-cylinder radius. 

Making further algebraic transformations on the equation © and taking into account equation 
d8]) we come to the following expression for the Goldreich- Julian charge density 



PGJ 2ttc Nrj 3 /(l) I 1 v 3 v)v 2 sm 2 6^' (U) 



Hereafter for the simplicity of calculations we assume that the inclination angle between the 
magnetic axis and the axis of rotation of the star is equal to zero. Figure CD presents the radial 
dependence of obtained general relativistic expression for the Goldreich- Julian charge density pcj 
normed in its Newtonian expression for several different values of NUT parameter. It can be found 
from the figure that even for comparatively small values of NUT parameter its influence on the 



-6- 



Goldreich- Julian charge density p G j CD} plays an important role. The value of the Goldreich- 
Julian charge density p GJ at the surface of the star is more sensitive to the NUT parameter. The 
tendency seen from the Figure is so that p GJ is vise proportional to the NUT parameter at the 
surface of the star and in asymptotics it reaches Newtonian expression. Here we took typical 
numbers for neutron star parameters as R = 10km, M = 2km and T = 0.1s. 

For the relativistic plasma the charge density p is proportional to magnetic field with the 
proportionality coefficient b eing constant along the given magnetic field line (see, for example, 
Muslimov & Tsvganl (119911) ) that is 

QB 1 f(r]) 



P 



(12) 



27TC NT] 3 f(l) 

where £ = 9/Q is the dimensionless angular variable, A(£) is an unknown function to be defined 
from the boundary conditions. One could insert expressions (fTTI) and (fT2l into the Poisson equation 
(OQ), and in the approximation of small angles 9 get the following differential equation 



R- 2 lN 



1 d 



-Air 



rf drj 
QB 1 







drj 

m 



2vrc AV f(l) 



1 



Nt] 2 9 



'if 



d_ 
89 



d_ 

89 



L I- 



7] J T]' 



l_a 2 1 

9d(f) 
+ A(0 



(13) 



Here we used a fact that in the small-angle limit 

4 sin 2 - 9 



sin 2 9 



1 . 



(14) 



Our further discourses are based on the extension of work iMuslimov & Tsvganl (|1992|) to 
NUT spacetime. Using dimensionless function F = 77$/$o> where $o = QB R 2 and variables rj 
and £, we can rewrite the equation ([T]) for the dimensionless electrostatic potential as 



dry 



ld_ ( d_ 



* (i - ? ) m 



1 - 



L I 



rj J r\" 



(15) 



where A(^) = [776(77) (1 - e/77) 1 / 2 ]" 1 . 

After performing Fourier-Bessel transformation 

and using relation 



i=l 



tF(r),{)MkitH 



E 



^ kiJi{ki) 



(16) 



(17) 
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one can obtain equation (1151) in the form 



1 - e - 

v 



/(I) 



kiJiyki 



L I 



(18) 



where jf = kfA 2 , ki are positive zeros of the functions J 



Considering a region near the surface of the star, where z — rj — 1 <C 1, and using following 
boundary conditions (that is the conditions of equipotentiality of the stellar surface and zero steady 
state electric field at r = R) 



i U=0 







dz 



2=0 



. 



(19) 



one can find the expression for the scalar potential $ near the surface of the star and corresponding 
to this poten tial component of the electr ic field En, being parallel to the magnetic field (see the 
discourses of Muslimov & Tsygan (| 1992 ) work): 



36$ r 



En 



£[K 



t=0 



exp 



36$ 
R 



i=0 



1 — exp 



kfJi(h) 



kj(l-y) 



kfJi(h) 



(20) 
(21) 



It should be noted that these formulae differ from the ones obtained in iMuslimov & Tsygan 



(119921) by that the parameter k gets replaced by (« — Le). 

Considering now the region O <C rj — 1 <C R LC /R, where \d 2 Fi/dr] 2 \ <C one can 



see that equation (1181) becomes 



rj 2 1 



/(I) 



kiJ\(ki 



K 



L I 



from which it immediately follows 



kf v '7(i) 



kiJ\[ki 



7,(1) 



77/77 



L 1- 



Ai 



(22) 



(23) 



Using this expression for Fj one can obtain the scalar potential in the region at distances 
greater than the polar cap size as 

V 



h el( K -Le) (i-l) (i-e 2 ) 



(24) 
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Corresponding to this potential component of electric field E» will look like 



11 Rdr] 



\£=constant 



FT A 2 



3(k — Le) 



a, 



(25) 



where E vac = (VlR/c)Bq is the characteristic N ewtonian value of the electric field generated near 
the surface of a neutron star rotating in vacuum iDeutschl (I1955|) . In Fig.|2]one can find the radial 
dependence of parallel electric field E\\ in terms of E vac for the different values of NUT parameter. 

The energy losses from the polar cap of the rotating star with n on-vanishing NUT-parameter 
can now be calculated. According to iMuslimov & Harding! (| 19971) the expression for the total 
power carried away by relativistically moving particles is 



L p = 2(-c J p$ dS) 



(26) 



In slow rotating NUT spacetime 



- p$ 



Air 



nB y R'eim 

N V 3 /(I) 



[( K _ L£ )(1_ K _ L (1_ £ ))](1_^ 



(27) 



Inserting (ITTT) into (l26l) and taking the integral one can get, that 

3 



(L 



p ) max 



[(k — Le) (1 — ac — L(l — e))\ E, 



rot ; 



where 



1 tt 4 B 2 R 6 



rot 



{E 



rot I Newt 



(28) 



(29) 



6 cV 2 (l) P{1 

and (E rot ) Newt is the standard Newtonian expression for the magneto -dipole losses in flat space 
time approximation. 



In the limiting case when / — > one could get the result IMuslimov & Harding! (| 19971) : 

3 



(L 



pjmax {1=0) 



-K(l - K)E. 



rot 



The ratio 



[L P ) 



p) max 



L(k + e-2ke) L 2 e(l-e) 



(L 



pjmax (i=0) 



k(1 — K.) 



k{1 — K] 



(30) 



(3D 



as function of NUT parameter is presented in Fig. [3j The dependence has a parabolic form. 
Namely, for small values of NUT parameter the energy losses are decreasing. Then the graph 
has a minimum and with further growth of NUT parameter, more energy is lost. Physically it is 
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due to the fact that contributions in accelerating electric field coming from the different parameters 
have the opposite signs. 

Taking into account that k = e{3 ~ e one can see, that the magnitude of additional terms 
arising from the NUT parameter is determined by the magnitude of L. For a millisecond pulsar 
with I ~ 10 3 cm, Q ~ 10s -1 and R ~ 10 6 cm we have 



d 



1 . 



(32) 



It shows that corrections to energy losses, concerned with non-zero NUT parameter cannot be 
neglected and in principle can provide an important information which will help to detect the 
gravitomagnetic charge. 

Equation (1281) has physical sense only if (L p ) max < E rot , that means the power of a polar cap 
accelerator in principle cannot exceed the total spin-down power of a pulsar if one suppose that 
electromagnetic radiation is powered by rotation. Using equations ( 1281 ), (|30b and (|3T| ) one can see 
that the physical sense is preserved for discussed magnetosphere model (when electric field created 
by rotation of star is dominating) if 



L 2 < 



3k(1 



1 . 



(33) 



Using the value k = 0.15 (|Muslimov & Tsyganl (119921) ) one can obtain upper limit for the value of 
NUT-parameter I < 1000cm. 



For practical useful applications the equation (1281) can be rewritten in terms of pulsar's ob- 
servable characteristics as the period P and its time derivative P = dP/dt: 



(PP) 

ma. 

where the expressions 
and 

have been taken into account. 



l[( K -Le)(l-K-L(l-e))]lj±p 



(PP) 



Newt 



(Lip)max 
(PP) 

Newt '- 



-/(fifi), 



2f_\ R 6 B 2 
3^ 



(34) 

(35) 
(36) 



In equation (1351) I is the general relativistic moment of inertia of the star (see e.g. lRezzolla & Ahmedov 



J d 3 Xv ^ye-* (r) pr 2 sin 2 9 



(37) 
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where e *W = l/y/—goo, p(r) is the total energy density, 7 is the determinant of the three metric 
and <i 3 x is the coordinate volume element. 

Period of pulsar and it's time derivative a re very prec i sely m easured quantities for a large 
number of pulsars (for example, in the paper of iKaspi et al.l (120061) t here is a P — P diagram for 



Arons 



(2007)). Thus, expression C34J for 



the 1403 catalogued rotation-powered pulsars, see also 
PP may indicate the possible existence and magnitude of NUT-parameter. The main difficulty 
encountered on this way nowadays is the uncertainty of estimation of the moment of inertia of the 
neutron star. But in future, when the moment of inertia of the neuron stars would be determined 
more precisely, one can in principle obtain the value of NUT-parameter from the observational 
data. 



Charged Particle Acceleration in a Polar Cap in Magnetosphere of Slowly Rotating 

NUT Star 



The origin of radio emission from the polar cap is one of the most mysterious, still unsolved 
question in the physics of pulsars. One of the likely scenarios is that particles are accelerated along 
open magnetic field lines and emit 7-rays that subsequently convert into electron-positron pairs 
under a strong magnetic field. The combination of the primary beam and pair plasma provides 
the radio emission mechanism. Due to this reason it is interesting to study particle acceleration 
conditions and equations of motion in the pulsar magnetosphere in the presence of NUT parameter. 



In our preceding paper iMofiz & Ahmedovl (120001) effects of general relativity on plasma 
modes along the open field lines of rotating magnetized neutron star have been studied. Here we 
investigate equations of motion of charged particle in the re gion just above the pola r cap surface 
of a neutron star with NUT parameter by extending results of lSakai & Shibatal (|2003l) to spacetime 
with gravitomagnetic charge. 



Equations of motion for particle with mass m, charge e, proper time r and 4-velocity 
dx^/dr look like 



m 



dr 



(38) 



where T^ x is the affine connection, = u^E v — u^E^ + rj^^UyBp is the electromagnetic field 
tensor, jf"^P is the Levi-Civita tensor. 



In the approximation of small angles 9 in the polar cap region the equations of motion ([38 
can be rewritten in the form 



N d 



s 2 ds V ds 



d<j>\ 1(1 + 1) 



D 



BqN 



J 

v~ 3 



(39) 
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(40) 



d _ 1 d(j) 

in the space-time outside the star described by Kerr-Taub-NU T metric in the s l ow ro tation limit 
©. It formally coincides with the equations of motion from ISakai & Shibatal (|2003r) . Here the 
Lorentz factor 7 = —v^u^ = Nv° = 1/(1 — V 2 ) 1 ^ 2 , the relative 3-velocity V 1 = v^h^/'y, the 
projection tensor h a!3 = g a/3 + ifu 13 and the following normalized variables 



2tiN{s)J{s) 

riB(s) 



const 



2ttB e 



mc* 



b(s) = — $(s) 
m 

27tN(s) Pgj (s) 
QB(s) 



(41) 



are introduced. To obtain equations (|39l) and ( |40l) the scalar potential $ was expanded in spherical 
harmonics Yi m (6, <fi) as $ = m <i(r)Y5 m (6>, <fi) and only the mode of the polar cap scale I ps 



7t/0 o was taken. In (1411) J is the electric charge current density, B = 
The approximate value of Goldreich- Julian charge density © is 



B| 



B r + O{0 2 



Pgj 



1 BQ 



1 



N 2 l 



(42) 



and consequently NUT-parameter will provide additional radial dependence to j. 



According to recent paper of lBeloborodovi (|2007h motion of charged particles in the vicinity 
os the polar cap of the pulsar is governed by equation (momentum of the particle is given in units 
of mc) 

dp 
dt 



eEn 



mc 



(43) 



And for parallel electric field E\\ one could get equation 

dEw 



dz 



4vr(p - p G j) 



pc ) 



(44) 



where r pc is the radii of the polar cap of the star. Rewriting (|44l) using d/dt = vd/dz, a = j/ c P gj, 



p = j/v and v = cp(l + p 2 



-1/2 



one could obtain (see iBeloborodovi (|2007l )) 



dE\\ 
~dT 



47TJ 1 



ap 



p A 



where 



a(z) = a - 



1 — k + Le 



I- k + Le){l + z/R)- 



(45) 



(46) 
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Then the system of equa tions for p(z) is graphically solved. Fig @] presents p(z) for several 
values of Le, when k — 0.15 (|Muslimov & Tsyganl (119921) ) and a = 0.999 and shows that influ- 
ence of NUT-parameter significantly changes the period of oscillations. From the Fig. [5] and Fig. [6] 
one could see when k — Le ~ 1 oscillations take place even for large a . For the values of k — Le 
very close to 1, the character of graphs almost has no dependence on the value of a . For old recy- 
cled neutron stars with almost zero angular momentum the effects connected with NUT-parameter 
may play single mechanism for pulsar radiation. 



Following lArons I (| 19981) and using equations (PTT)) . (fT2j) and (IT4|) . one can obtain for the 
accelerating potential drop: 



while in Newtonian case it looks like 



(k - Le)<& 



pole 



pole 



1- ^ 



— ) 

PbJ 



(47) 



(48) 



where ps is the radius of curvature of the magnetic field lines, $ po z e 
magnetic moment. 



Vt 2 n/c 2 and fi is the 



Taking into account the frame dragging effect noticeably improve comparison between theory 
and observations. As R/ pB ~ 10~ 2 P -1 / 2 and k ~ 10P" 1/2 , the energy of curvature gamma 
ra ys in relativistic frame - dragg ing case may rise up to be 1000 times greater than that occurs 



in 



Arons & Scharlemannl (|1979|) pair creation theory (see lArons I (|1998|) ). As the effect of NUT- 



parameter seems to have the same order as the effect of frame dragging, it makes an additional 
contribution to the energy of curvature gamma photons and thus may play noticeable role in the 
formation of the plasma magnetosphere of neutron star with gravitomagnetic charge. 



4. Conclusion 

We have considered astrophysical processes in the polar cap of pulsar magnetosphere in as- 
sumption of slowly rotating NUT space-time. In particular, general-relativistic corrections to the 
Goldreich- Julian charge density, electrostatic scalar potential and accelerating component of elec- 
tric field being parallel to magnetic field lines in the polar cap region due to the presence of gravit- 
omagnetic charge are found. The presence of NUT-parameter slightly modulates Goldreich- Julian 
charge density near the surface of the star. However as it is known the effective electric charge den- 
sity i.e. difference between Goldreich-Julian charge density pcj (being proportional in the case of 
flat space-time to Ct ■ B) and electric charge density (being proportional to B) in the star magneto- 
sphere is responsible for the generation of electric field being parallel to magnetic field lines. This 
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difference is equal to zero at the surface of the star and changes with the distance from it due to the 
fact that p can not compensate p GJ . General relativitistic terms arising from the dragging of inertial 
frames and presence of the gravitomagnetic charge give very important additional contribution to 
this difference. Both of these terms depend on the radial distance from the star as 1 /r 3 and have 
equally important influence on the value of accelerating electric field component generated in the 
magnetosphere near the surface of the neutron star. 

These results are applied to find an expression for electromagnetic energy losses along the 
open magnetic field lines of the slowly rotating NUT star. It is found that in the case of non- 
vanishing NUT-parameter an additional important term to the coefficient of standard magneto- 
dipole energy losses expression appears. Comparison of effect of NUT-parameter with the already 
known effects shows that it can not be neglected. Obtained new dependence may be combined with 
astrophysical data on pulsar periods slow down and be useful in further investigations on possible 
detection of the gravitomagnetic monopole. 

It is also shown that the presence of gravitomagnetic charge has influence on the conditions of 
particles motion in the polar cap region. From derived results it can be seen, that NUT-parameter 
modulates the period of oscillations of particle's momentum. As the effect, connected with NUT- 
parameter has the angular velocity of the star's rotation in denominator, it will be increased for 
the old stars, where the frame-dragging effect will be decreased. Thus possible candidates for 
verification of the existence of NUT parameter are old compact objects with extremely low rotation 
period. 
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Fig. 1. — Radial dependence of the Goldreich- Julian charge density normed in its Newtonian 
expression for the different values of the NUT parameter. 
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Dependence of Em / E vac on NUT parameter 
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Fig. 2. — Radial dependence of the accelerating component of electric field normed in its vacuum 
magnitude for the different values of the NUT parameter. 
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Dependence of (L p ) m ax/C-p)max, (i=o) on N UT parameter I 
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Fig. 3. — Ratio of polar cap energy losses into its zero gravitomagnetic charge expression as a 
function of NUT parameter. 
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Fig. 4. — The dependence of the momentum of charged particle, extracted from the polar cap by 
E\\, on height for different values of NUT-parameter. The range of used values: B = 3 x 10 12 G, 
P — Is, k — 0.15, a = 0.999, rotator is assumed to be aligned. Considered region is z < 100m 
which lies within the polar cap radii. 
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Fig. 5. — The dependence of the momentum of charged particle, extracted from the polar cap by 
E\\, on height for a — 1.1, k — and different values of NUT-parameter. 
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Fig. 6. — The dependence of the momentum of charged particle, extracted from the polar cap by 
E\\, on height for a = 2, k = and different values of NUT-parameter. When NUT-parameter Le 
is very close to 1, oscillations take place independently on the value of a . 



